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Acronyms and Abbreviations 


aluminum 

software application used for spacecraft MMOD risk assessments 
Ballistic Limit Equation 

Earth Entry Vehicle 

Element Identification Number (unique identification number) 
Engineering Services Contract Group 

Finite Element Model 

Hypervelocity Impact Technology Facility 

Hypervelocity Impact 

Hypervelocity Impact Technology Group 

Integrated Design, Engineering, Analysis, and Simulation software 
International Space Station 

Johnson Space Center 

Meteoroid Environment Model 

Multi-Layer Insulation 

Micro-meteoroid Orbital Debris 
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National Aeronautics and Space Administration 

Orbital Debris Environment Model 2000 

Property Identification Number (used to index an element to a specific BLE) 
Probability of No Penetration 
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Thermal Protection System 


White Sands Test Facility 


1. Introduction 


This report provides results of a Micrometeoroid 
and Orbital Debris (MMOD) risk assessment of 
the Mars Sample Return Earth Entry Vehicle 
(MSR EEV). The assessment was performed using 
standard risk assessment methodology illustrated 
in Figure 1-1. Central to the process is the Bumper 
1-2), which 
calculates the critical penetration risk based on 
geometry, and flight 
parameters. 


risk assessment code (Figure 


shielding configurations 
The assessment process begins by building a finite 


element model (FEM) of the spacecraft, which 
defines the size and shape of the spacecraft as well 
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Figure 1-1. MMOD Risk Assessment Flowchart 


as the locations of the various shielding configurations. This model is built using the NX I-deas software package from 
Siemens PLM Software. The FEM is constructed using triangular and quadrilateral elements that define the outer shell of the 
spacecraft. Bumper-II uses the model file to determine the geometry of the spacecraft for the analysis. 


The next step of the process is to identify the ballistic limit characteristics for the various shield types. These ballistic limits 
define the critical size particle that will penetrate a shield at a given impact angle and impact velocity. When the finite 
element model is built, each individual element is assigned a property identifier (PID) to act as an index for its shielding 
properties. Using the ballistic limit equations (BLEs) built into the Bumper-II code, the shield characteristics are defined for 


each and every PID in the model. 


Spacecraft Configuration (I-DEAS Finite Element Model) 
¢ Describes spatial relationships of spacecraft components 
¢ Defines spacecraft orientation (velocity and zenith directions) 


¢ Defines M/OD shield regions ¢ Shadowing 


¢ 90 debris threat cases 2 
and 149 meteoroid threat cases 


¢ Approximately 120,000 elements in ISS assembly complete mated configuration FEM 


Critical Particle Diameter Calculation (RESPONSE) 


¢ Protection capability 


assessed for each element in the FEM 


Meteoroid & Debris Environments (GEOMETRY) 
¢ Threat directions aah 
¢ Velocity distribution 


Threat Elements 


__. Meteorord 


meteoroid and debris 


Computation of Penetrating Flux and PNP (SHIELD) 
Graphical Interpretation of Results (EXCEL & I-DEAS) 


Space Station Orbital Debris Threat Assessment 
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Figure 1-2. BUMPER-IIT MMOD Risk Assessment Code 


impact Aik From 1mm Debris 


0.985074 
0.985522 
0.997443 


0.934622 


Debris Penetration Risk 





The final stage of the analysis is 
to determine the probability of 
no penetration (PNP) on the 
spacecraft. This is done using 
the micrometeoroid and orbital 
debris environment definitions 
that are built into the Bumper-II 
code. These engineering models 
take orbit 
inclination, altitude, attitude and 
analysis date in order to predict 
an impacting particle flux on the 


into account 


spacecraft. Using the geometry 
and shielding characteristics 
previously defined for the 
spacecraft and combining that 
with the 


environment model calculations, 


information 


the Bumper-II code calculates a 
probability of no penetration for 
the spacecraft. 


This report is prepared for NASA Ames Research Center (POC: Ethiraj Venkatapathy), and responds to a statement of work 
that is provided in Appendix A - MSR EEV MMOD Impact Risk Analysis Statement of Work. 
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Z. Requirements 


The results of this analysis were compared with a probability of containment not assured due to MMOD impact of 4.0x107’. 
This value was assigned during the August 9, 2011, monthly status teleconference. Project engineers made the allocation 
from the probabilistic risk assessment (PRA) fault tree (Gershman et al. 2005), reproduced below as Figure 2-1. MMOD 
impacts of the EEV can only result in a fault through the “EEV fails during correctly targeted entry” branch of the tree. 
Therefore, an allocation 4.0x10-’ to MMOD impact leaves nothing for other modes of failure during entry. Thus, unless a 
future reliability analysis shows that there are no other modes of failure occurring in this branch, the 4.0x10~’ probability 
should be taken as an upper bound on the requirement. 


CONTAINMENT P=13E6 
NOT ASSURED 





EXTERNAL EEV EEV INCORRECTLY FEV FAILS DURING LOSS OF ERV IMPACTS 
CONTAMINATION TARGETED CORRECTLY CONTAINMENT EARTH 
TARGETTED DURING OR AFTER CONTAMINATED OR 
ENTRY LANDING WITH EEV 





ATTACHED 
4E-7 3E-7 1E-9 


1E-9 6E-f 


Figure 2-1. MSR EEV - PRA Fault Tree. 
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3. MSR EEV Thermal Protection System Design 


The May 10, 2000 design for the MSR EEV is shown in cross section in Figure 3-1. The EEV is designed to land without 
parachutes, for maximum reliability. The blunt shape slows the EEV to a terminal speed of 41 m/s (92 mph) and the Earth 
impact loads are taken up by the spherical crushable structure. The orbital sample container is visible Figure 3-1 at the center 
the energy-absorbing sphere, which is composed of carbon foam plus webs. The blunt heat shield is a thin carbon-carbon 
Structure; filled with carbon foam at the center and void at the circumference. 


The current MSR EEV thermal protection system (TPS) has two parts. The forward TPS consists of 12mm of fully dense 
carbon phenolic. Carbon phenolic was selected for high reliability because of the extensive flight heritage and constitutes 1/3 
of the total mass of the EEV. The aft TPS material has not been selected yet. Weight estimates were made with 10 mm of 
SLA-561V, and that is what is used in the MMOD impact risk assessment. Multilayer thermal insulation covers the MSR 
EEV’s forward TPS and a portion of the aft TPS. 
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Figure 3-1. MSR EEV Shield Configuration. 
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4. Analysis 


This section provides details on the analysis procedure and the Bumper-II input parameters used in the MSR EEV analysis. 


4.1 Failure Criteria 


Failure was defined as a perforation of the TPS down to the underlying structure. This failure criterion corresponds to a 
perforation through the 10 mm SLA-561V aft TPS or a perforation through the 12 mm carbon phenolic forward TPS. 


These criteria are less conservative than what are used to evaluate the Orion MPCV. Penetration of no more than 25% to 80% 
of the Orion heat shield TPS thickness is allowable, depending on the location, and no more than 62% to 90% of the back 
shell TPS (Bohl, 2011). If the more conservative Orion TPS failure criteria apply to EEV, then the risk of containment not 
assured was underestimated in this report. Further work is recommended to confirm or change the EEV TPS failure criteria. 


4.2 Configuration 


Robert Maddox provided an EEV CAD model to JSC/KX engineers on April 14, 2011. A rendering of the CAD model is 
shown in Figure 4-1. 





Figure 4-1. MSR EEV CAD Model. 


Kevin Hoffman created the shell finite element models from the EEV CAD model for the MMOD analyses. Two components 
of the EEV were modeled: the forward and aft TPS systems. These two systems were divided into a total of 13 distinctive 
regions based on the structure below the TPS. Figure 4-2 is a rendering of the EEV FEM with the 13 regions colored and 
labeled. Dividing the model into regions permits the assignment of unique failure criteria to each region. Should further test 
and analysis show that penetration down into the structure is allowable, different failure criteria could be assigned to the 13 
regions. 





Aft Side Forward Side 


Figure 4-2. MSR EEV Analysis — Finite Element Model. 
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However, for this analysis each of the 13 regions was given one of two failure criteria: penetration through the thickness of 
the forward TPS or penetration through the thickness of the aft TPS structure. Furthermore, should it be necessary to use 
failure criteria similar to the Orion TPS, then the failure criteria will vary from the windward to the leeward sides, and it will 
be necessary to add more regions to the finite element model. 


The I-DEAS finite element model variable that is mapped into the failure criterion used by BUMPER-II is called a property 
ID (PID). Thus, each of these 13 risk regions called a PID in this report. Some of the geometric characteristics of the 13 PIDs 
are listed in Table 4-1. The first column lists the PID number used in the I-DEAS universal file. The second column lists the 
label assign to the PID by the finite element modeler. The third column lists whether the PID was assigned an aft TPS failure 
criterion or a forward TPS failure criterion. The fourth column lists the color of the PID in Figure 4-2. The fifth, sixth and 
seventh columns list the starting element ID for the PID, the ending element ID, and the total number of elements in the PID, 
respectively. The eighth column lists the surface area of the PID in square meters. 


Table 4-1. MSR EEV Analysis - Risk Region, PID, EID and Area Map. 
MMISR EEV FEM - Risk Region, PID, EID, and Area Map 


Element ID Range Surface 
Risk Region TPS Side Color Sd eal aa Area (m’) 


Shadowing white 
LID TPS (No-MLI) Aft orange 
Aft TPS (Body Foam, No-MLI) Aft red 
Aft TPS (Gap, No-MLI) Aft pink 
Aft TPS (MLI, Gap) Aft blue 
Fwd TPS (MLI, Gap) Forward yellow 
Fwd TPS (MLI, Body Foam) Forward cyan 
Fwd TPS (Double MLI, Gap) Forward | magenta 
LID Push Pad Aft It blue 
LID Pin Aft yellow 
Body Mount Bolt Aft green 
Unknown 1 Aft white 
Unknown 2 Aft white 
Unknown 3 Aft white 


The EEV is transported to Mars and then back to the vicinity of Earth by a carrier vehicle. This is illustrated in the graphic 
shown in Figure 4-3. The 2018 launch on the left hand side of the plot lands a rover on Mars to collect samples and 
rendezvous with the Mars ascent vehicle. The 2022 launch sends the MSR orbiter to Mars with the EEV. The third and final 
launch in 2024 sends the MSR lander with ascent vehicle to Mars. The lander inserts a specimen (the orbital sample or OS) in 
the ascent vehicle and the ascent vehicle rendezvous with the orbiter. The orbiter then leaves Mars and releases the EEV in 
the vicinity of Earth for Earth entry. 
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Figure 4-3. MSR EEV Analysis - Mission Configuration (Mattingly and May, 2011). 


This orbiter is much larger than the EEV as shown in the upper right inset of Figure 4-3 with the dark yellow EEV mated to 
the side of the orbiter. The large orbiter shadows the EEV from the meteoroid environment for the majority of the mission, 
reducing the risk of meteoroid impact on the aft TPS. A notional model of the orbiter was constructed using a cylinder as 
shown in Figure 4-4. 





Figure 4-4. MSR EEV Analysis - KEV FEM with Transit Vehicle. 
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Additionally, a FEM with the EEV completely enclosed in an MMOD “garage” (for all mission phases up to Earth reentry) 
was used for a concept study and is shown in Figure 4-5. 


MMOD Shield concept has a 
35 cm standoff from the EEV 


Transparent view 





Figure 4-5. MSR EEV Analysis - MMOD Shield Concept FEM. 


4.3, Assessment Inputs 


MSR EEV analysis parameters are summarized in Table 4-2. Version “CEV 1.73” of the BUMPER-II code was used for the 
assessment. The current orbital debris engineering environment, ORDEM2000 was used with a constant projectile density of 
2.8 grams per cubic centimeter. The MEM environment was used for meteoroid risk calculations. [IP MEM v2.1 was used for 
the transit phases and MarsMEM was used for the orbital phase. Additional information is available in section 4.4.3. 


Table 4-2. MSR EEV Analysis Parameters. 
MSR EEV Analysis Parameters 


Analysis Code BUMPER-II CEV v1.73 


Bea easlhLa 


Debris Density 


steseneasenceceececeecessereceecececsacescecesescescecesea Qeeseasessesessraeseacesceseaceatesnaee: 


Number of Debris Threat Cases 


Number of Meteoroid Threat Cases 





The MSR EEV mission consists of 15 phases (Table 4-3) from launch of the MSR EEV (with transit vehicle) to reentry of the 
EEV only. The total exposure duration of the MSR EEV to the MMOD environments is approximately 4.83 years. For this 
analysis, meteoroid impacts during the launch/LEO and EEV reentry phases were ignored. Additionally, only orbital debris 
impacts during EEV reentry were considered. 
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Table 4-3. MSR EEV Analysis - Mission Phases. 


MSR EEV Mission Phases 
1 Launch-LEO 0.00023 2 hrs, ignored meteoroids and OD 
Earth to Mars transit 0.70833 | 
3-11 | Mars Arrival (aerobraking) 0.41096 











| 12 | Mars Orbit 3.00000 
Mars departure orbit 0.00548 
| 14 | Mars to Earth transit 0.70833 
15 | EEV Entry 0.00023 2 hrs, assessed OD, ignored meteoroids 


pj total 88856 | 


Due to the preliminary status of the MSR EEV mission design, the attitude of the MSR transit vehicle has not been 
determined. Therefore the approach taken for the MSR EEV analysis was to use 24 attitudes (90° rotations about the axes of 
the spacecraft coordinate system) for the mission phases while the EEV is attached to the transit vehicle. These 24 attitude 
cases were used only with the meteoroid environment (orbital debris only considered during EEV Earth reentry) and are 
summarized in Table 4-4. The rotations are expressed in a “Bumper Roll, Pitch, Yaw” or “RPY” space fixed sequence, where 
“Roll” is a rotation about the “X” axis of the reference frame, “Pitch” is a rotation about the “Y” axis and “Yaw” 1s a rotation 
about the “Z”. The XYZ reference frame for the attitude cases is illustrated in Figure 4-6. 











Table 4-4. MSR EEV Analysis — Meteoroid Attitude Cases. 


MSR EEV Attitude Cases 


coset ae 
case# [xX |  Y__ 


© 
OO! 0O?0:0:!:0:!:0:0O:O 


0 


O20! 0: 0! 01010! 0O!10O:0:0:0:0 





Additionally, the 24 attitude cases are depicted graphically in Figure 4-6. Further information regarding the coordinate 
system used for the meteoroid environment can be found in Appendix E - MEM Coordinate System. 
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Transfer Orbit 
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Place picture from left in 
dotted square above 


Figure 4-6. MSR EEV Analysis — 24 Meteoroid Attitude Cases. 


For the EEV Earth reentry phase (orbital debris only), a single attitude (Figure 4-7) was considered with the EEV pitched 
down 20° with respect to the local horizontal. The selection of the Earth reentry trajectory and attitude is discussed in 


Section 4.4.1.5. 
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Local Horizontal 


20° pitched 
down wrt local 


horizontal 
EEV velocity vector 





Figure 4-7. MSR EEV Analysis - EEV Earth Reentry Attitude. 


4.4 Environment Models 


4.4.1 Trajectory 


The spacecraft trajectory through the orbital debris and the meteoroid environments is one of the factors that determine the 
penetrating fluence of particles. The MSFC Meteoroid Environment Office used recent Mars spacecraft trajectories to bound 
the range of possible meteoroid fluence. Bill Cooke’s discussion of the reasoning used to select trajectories for the meteoroid 
impact analyses is reproduced below from a June 2011 email. 


The time of flight in low Earth orbit for both outbound and inbound trips is approximately 1.75 hours, which is not a 
significant amount of time exposed to the meteoroid environment. It is not necessary to model these phases of the mission for 
the number of impacts by meteoroids analysis. (However, it is necessary to analyze these phases to determine the orbital 
debris fluence. While the orbital debris fluence is small, the mean number of impacts requirement is orders of magnitude 
smaller than the requirement for other spacecraft analyzed by JSC/KX engineers to date. The discussion of the earth entry 
trajectory used in the number of impacts by orbital debris particles appears in Section 4.4.1.5.) 


The exact trajectory for the Earth-Mars-Earth interplanetary cruise is unknown. Several NASA Mars missions were 
investigated and their interplanetary flight paths were used to generate meteoroid environment files and bound the possible 
variations in flux due to unknown flight path. The following mission trajectories were readily available from the JPL 
Horizons website in the necessary format: Mars Exploration Rover — Spirit, Mars Exploration Rover — Opportunity, Mars 
Reconnaissance Orbiter (MRO), and Mars Odyssey. Mars Pathfinder orbital elements for half of the interplanetary cruise 
were located on the mission webpage and used for comparison. Also, a simple Earth-Mars-Earth Hohman trajectory was 
generated for comparison. The Spirit and Opportunity missions represent “type 1” interplanetary transfers, the fastest 
possible. The cruise phases lasted 5 months. MRO and Odyssey had longer cruise phases of 7 months and a Hohman 
minimum-energy transfer represents the longest case with a cruise time of 8.5 months. 


Since the meteoroid environment directionality is fixed with respect to the Sun, the spacecraft can be exposed to different 
sources of meteoroids depending on its flight path and attitude. It is assumed this spacecraft is not rotating and the ram 
surface (X) is defined by the spacecraft’s velocity vector, port surface (Y) is defined by the angular momentum vector and 
zenith surface or space facing surface (Z) is the cross product of the velocity vector and angular momentum vectors. 


Upon insertion into Mars orbit, the spacecraft spends 150 days aerobraking into a circular orbit. This aerobraking phase was 
simulated by generating orbits at various apoapse altitudes but with the same inclination, node and periapse altitude. It is not 
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known how long the spacecraft will spend at each altitude so an average environment was constructed over one orbit at each 
altitude during this aerobraking phase. 


The spacecraft trajectories from JPL Horizons used for the Earth to Mars and the Mars to Earth transit phases and the 
Keplerian elements used for the elliptical orbit phases are listed in Table 4-5. These choices are described in greater detail in 
what follows. The first column of Table 4-5 lists the mission phase ID number, the second the mission phase label, the third 
the trajectory type, the fourth the Mars apoapse for mission phases 3 to 13 and the Earth apoapse for mission phase 15, the 
fifth column the Mars periapse for mission phases 3 to 13 and the Earth periapse for mission phase 15, the sixth column lists 
the right ascension of ascending node, the seventh column lists the argument of perigee, the eighth column the orbital period 
and the ninth column the inclination of the orbit. 


Table 4-5. Spacecraft Trajectories and Keplerian Elements for the 15 Mission Phases. 
Mission phase Trajectory Apoapse’ Periapse Node 0 Period | Inclination 
Cy Cry (deg) __ (hrs) (deg) 


2 Earth to Mars transit Spirit, MRO, Hohman 


4 = Mars Arrival (aerobraking) 1 orbit 49,811 3,897 152 37 45 
10 Mars Arrival (aerobraking) 1 orbit 4,087 3,897 152 45 


15 EEV Entry Genesis (EI-2 hrs to EI) 1,382,034 — 6,378 ae 90 1,599 oe 





4.4.1.1 Launch Phase 


An MSR reference trajectory for launch is unavailable for use in generating an average environment file. However, by 
definition, a launch trajectory is a power track, and hence hyperbolic in Earth centered inertial (ECI) coordinates. This was 
confirmed by using JPL Horizons to calculate the MRO launch trajectory. The launch phase covers the first 2 hours of 
mission elapsed time. 


4.4.1.2. From Earth to Mars Cruise Phase 


An MSR reference trajectory to Mars is unavailable for use in generating an average environment file; however several 
recent NASA missions have their interplanetary cruise trajectories available on the JPL Horizons website. In order to 
describe and bound the meteoroid environment for this phase of the mission, three trajectories were used to develop the 
meteoroid flux-direction files. The Mars Exploration Rover — Spirit trajectory was chosen because it followed a “type 1” 
interplanetary transfer, the fastest possible, arriving at Mars in 5 months. The MRO spacecraft trajectory was chosen because 
it had a 7 month interplanetary cruise approaching Mars from the south and was typical of an average transit time. Finally, a 
minimum-energy Hohman transfer trajectory was used where the cruise phase lasted about 8.5 months. 


4.4,1.3_ Aerobraking, Mars Orbit and Mars Departure Phases 


While some information was available for the initial orbit, other information had to be assumed in order to approximate a set 
of environment files for several different orbit orientations. There are 11 meteoroid environment files that describe the 
meteoroid environment during the aerobraking phase, Mars orbit and Mars departure orbit phases. The aerobraking phase 
occurs over 150 days, 10 orbits were constructed to simulate the falling apoapse. The last file represents the final orbit where 
the vehicle is presumed to spend the most time before raising the apoapse for the Mars departure orbit. Each file represents 
the “average” and | sigma environment over | orbit at that altitude sampled in 60 second time steps. 


The orbit parameters at Mars were described by the following: an initial orbit with period = 48 hours, periapse altitude = 500 


km, inclination = 45 degrees, Node-IAU = 15.6 degrees, argument of periapse = -0.2 degrees, and Mean anomaly = 0 degrees 
at 12:00 ET on 7/10/2023. IAU is a vector defined in a Mars centered coordinate frame. The [AU- vector is defined to be 
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along the intersection of the Mars Mean Equator of Date plane and the Earth Mean Equator of Epoch J2000 plane, and is 
positive in the direction of the ascending node of the Mars mean equator of date on the Earth mean equator of epoch J2000. 
This vector can be found using a relationship between the vernal equinox vector and the angle between the IAU vector and 
equinox. This angle is about 136 degrees. Substituting the IAU angle in the above equation we can solve for the Node 
location of the orbit plane which is around 152 degrees. This IAU vector and Mars coordinate frames are described in an 
interoffice memorandum 343B-2006-004 dated August 15, 2006 from NASA JPL. An STK scenario was developed using 
Astrogrator’s ““Targeter” to simulate the initial Mars orbit, aerobraking phase, and final orbit mission for MSREEV. State 
vectors were found in a Mars centered inertial J2000 reference system for input into MarsMEM, for each orbit segment 
described above. Once the aerobraking phase ends the spacecraft resides in a more circular orbit with altitude of 572 km and 
inclination of 45 degrees, the node is “TBD”. The orbiter then rendezvous with the Mars surface sample at an unknown 
altitude between the ranges of 522-572 km for both periapse and apoapse distances. Sometime over the next 3 years between 
December 2023 and July 2026, the orbit period is raised to somewhere between 2 hours and 48 hours and then over a period 
of a few days in July 2026, the periapse is lowered to 300 km and apoapse is raised (if not already there) to a period of 48 
hours just before the trans-Earth insertion maneuver. 


4.4.1.4 From Mars to Earth Cruise Phases 


For the interplanetary return phase of the mission, there were no specific trajectories described, and no recent NASA Mars 
missions to model this. A Hohman transfer orbit was again used to describe the state vectors over the return trip. States were 
sampled in | day time increments over 8.5 months. 


4.4.1.5 EEVEntry 


The Genesis trajectory was selected as an example of an Earth return vehicle reentry trajectory. The JPH Horizons web site 
was used to generate the Keplerian elements in Earth centered coordinates of the last two hours of the Genesis trajectory. The 
header for the Horizons output file documents the input parameters used to generated the trajectory and is included as Section 
9 of this report. The Keplerian elements from Horizons are plotted in Figure 4-8. The elements are constant until the last 8 
minutes before entering the Earth’s atmosphere (Earth interface), so the constant values were used for the whole trajectory. 
These are listed in Table 4-5. (Two hours corresponds to 7,200 seconds, during which time the Genesis sample return capsule 
traveled from 34,000 km altitude to Earth interface. Geosynchronous orbit is about 36,000 km altitude, so for the last two 
hours of the trajectory the Genesis Earth return capsule was within the Earth orbital debris belt.) 
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Figure 4-8. Keplerian elements from JPL Horizons for the last 2 hours of the Genesis orbit before Earth Interface. 


The spacecraft attitude during reentry was estimated by assuming the apex of the forward TPS cone was pointed in the 
spacecraft velocity direction. The scalar product of the spacecraft position vector and the spacecraft velocity vector was then 
computed to obtain the cosine of the angle between the position vector and the velocity vector. The spacecraft position vector 
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is normal to the local horizontal, hence the angle between the spacecraft position vector and the spacecraft velocity is the 
complement of the angle the spacecraft velocity makes to the local horizontal. This angle is cross plotted with the spatial 
density of | mm orbital debris particles during the year 2005 for altitudes above the equator. The maximum spatial density 
occurs for angles with respect to the local horizontal around 20 degrees, which was used for the attitude over the whole of the 
Earth entry phase. 
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Figure 4-9. EEV attitude as a function of altitude cross plotted with the spatial density of 1 mm diameter orbital 
debris particles as a function of altitude. 


4.4.2 Orbital Debris Environment 


The ORDEM2000 orbital debris environment model (Liou, et al. 2002) was used for the Earth entry phase of the analysis. 
ORDEM2000 is the latest released version from the JSC Orbital Debris Program Office (ODPO) and is valid up to 2000 km 
altitude. Bjorkman (with Matney’s concurrence) extended the ORDEM2000 environment to larger altitudes by extrapolating 
the orbital debris spatial density at 2,000 km using an inverse square dependence on radial distance from the center of the 
Earth. This extrapolation is apparent as the smooth function between 2,000 and 36,000 km in Figure 4-9. The extrapolation 
between 2,000 and 36,000 km altitude is thought conservative, perhaps by a wide margin. The degree of conservatism that 
the extrapolation has added to the results is unknown. The soon to released ORDEM 3.0 environment will calculate flux on 
spacecraft from LEO to GSO, and could provide a check on the analyses performed for the EEV. 


The orbital debris flux varies with the solar activity during the 11 year solar cycle. Years of high solar activity heat the 
atmosphere and increase the rate of orbital debris orbital decay. These years have a small flux. Years of low solar activity 
have less solar heating of the atmosphere, less orbital decay, and relatively larger fluxes. The oscillation of the orbital debris 
flux with the 11 year solar cycle is superimposed on a more uniform rate of flux increase resulting from the accumulation of 
material in Earth orbit from the historically more-or-less constant launch rate. Since the MSR Earth return vehicle is assumed 
to reenter Earth’s atmosphere in the year 2026, these two factors involve a 25 years extrapolation from the creation date of 
the ORDEM2000 environment model. The effect of changes in the period of the solar cycle and the historical launch rate on 
the ORDEM2000 prediction of the orbital debris flux in 2026 have not been considered. 


A third, and more restrictive, ORDEM2000 use limitation is that it only calculates the orbital debris flux on spacecraft in 
closed elliptical orbits around the Earth. The MRO launch trajectory was computed in ECI coordinates with JPL-Horizons, 
but was found to be hyperbolic, hence not suitable for use with ORDEM2000. The Genesis sample return capsule entry 
trajectory was calculated with JPL-Horizons and was found to be elliptic and could be used with ORDEM2000. Thus this 
EEV analysis only considers orbital debris impacts during the Earth entry phase and not during launch. This may not be a 
significant oversight if meteoroid shields are added to the EEV bio-shield. If they are, then the risk of EEV TPS failure due to 
orbital debris impact may be significantly smaller than the risk during entry, when there is no shield. If it isn’t possible to add 
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meteoroid shields to the bio-shield, then the TPS risk of failure during launch, may be as large as during entry, and shouldn’t 
be overlooked. 


Bjorkman recently modified ORDEM2000 to output the flux for segments of an elliptical orbit. This modification was used 
to calculate the flux during the last 2 hours of the Genesis entry using the Keplerian elements listed in Table 4-5. However, 
ORDEM2000 is coded to divide an elliptical orbit into no more than 9,999 segments, so the last two hours of the extremely 
elliptical Genesis entry orbit covered only 13 segments. The cross sectional area flux for those 13 segments were written out 
and cross plotted with the spatial density of 1 mm size orbital debris particles during the year 2005 above the equator. 
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Figure 4-10 . Flux of 1 mm OD particles on the EEV as a function of altitude cross plotted with the spatial density of 1 
mm diameter orbital debris particles as a function of altitude. 


The plot illustrates two issues. First, 12 of the 13 points lie at altitudes above 2,000 km, and these 12 points add up to 68% of 
the orbital debris fluence on the EEV during entry. The flux from altitudes above 2,000 km was computed from an 
extrapolation of the spatial density. The extrapolation function was selected because it was thought to give a conservative 
answer. The degree of conservatism is unknown and won’t be quantifiable until the ORDEM2010 environment is released. 
The second issue is the spatial density rapidly varies with altitude below 2,000 km and the flux for altitudes below 2,000 km 
has been approximated by only one flux; the flux at 820 km altitude. 


Regardless of the issues, estimates of the impacting and penetrating flux were made with the existing ORDEM2000 orbital 
debris environment model. Should the ORDEM2000 results raise issues with the EEV design, then it may be prudent to 
modify the ORDEM 3.0 orbital debris environment model for EEV and recalculate the probability of TPS failure before 
modifying the EEV design. 


4.4.3. Meteoroid environment 


The MSFC Meteoroid Environment Office (MEO) provided a description and 13 models of the meteoroid environments for 
the applicable mission phases (Cooke, 2011). (See (McNamara et al, 2004) for a description of the MEM computer program 
used to generate the environments.) The mission phase number, mission phase label, meteoroid environment file name, and 
cumulative cross sectional area flux of 1 microgram or larger meteoroids are listed in columns 1 to 4 of Table 4-6, 
respectively. 


Table 4-6. Environment Files Used for the Fifteen Mission Phases. 
MEM.DAT or Cross sectional area flux of 1 microgram 
Mission phase ORDEM2K.DAT meteoroids 
file name (#/sq m/yr) 


Launch-LEO N/A 
Earth to Mars transit AvgOutIPMEMIgloo.out 4.138+0.570 
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MEM.DAT or Cross sectional area flux of 1 microgram 
Mission phase ORDEM2K.DAT meteoroids 
file name (#/sq m/yr) 


Mars Arrival 
(aerobraking) 
Mars Arrival 
(aerobraking) 
Mars Arrival 
(aerobraking) 
Mars Atrival 
(aerobraking) 
Mars Arrival 
(aerobraking) 
Mars Arrival 
(aerobraking) 
Mars Arrival 
(aerobraking) 
Mars Arrival 
(aerobraking) 


ace abo9_AvgMEMIelooDist.out 1.45840.247 


12 Mars Orbit abol0_AvgMEMIglooDist.out 1.458+0.246 
13 Mars departure orbit FinalOrbitAvgMEMIglooDist.out 1.662+0.300 
14 Mars to Earth transit AvgInIPMEMIgloo.out 3.898+0.603 


15 EEV Entry ORDEM?2K-rentry.dat 


abol_AvgMEMIglooDist.out 1.656+0.278 
abo2_AvgMEMIglooDist.out 1.658 0.278 
abo3_AvgMEMIglooDist.out 1.657+0.277 
abo4_AvgMEMIglooDist.out 1.656+0.276 
aboS_AvgMEMIglooDist.out 1.649+0.274 
abo6_AvgMEMIglooDist.out 1.646 0.274 
abo7_AvgMEMIglooDist.out 1.630+0.271 


abo8_AvgMEMIglooDist.out 1.451+0.245 





The cumulative cross sectional area flux values listed in column 4 as a function of mission phase reflect the decrease in 
spatial density of meteoroids with heliocentric radial distance. The cumulative cross-sectional flux of meteoroids in Earth 
space at 1 AU is approximately 5.95 m~yr'' for particles down to 1 microgram. At Mars this flux drops to about 2.34 m7yr’ 
with an average speed of 18.2 km/s. 


The cross sectional area flux listed in Table 4-6 1s listed with the estimate of standard deviation made by engineers at the 
Meteoroid Environment office. The method used to estimate the standard deviation is discussed in Section 4.4.3.4. 


The environment files are required because the background or “sporadic” meteoroid environment is not isotropic. Each file is 
a table of flux values for a distribution of approach angles. The distribution of approach angles comes from the 6 sources of 
meteoroids. Comets contribute the majority of the meteoric material in the inner solar system and asteroids contribute a much 
smaller portion. Within 1.5 AU, the sporadic sources appear fixed in their direction about the Sun. In the heliocentric ecliptic 
frame, the most active sources are the helion and anti-helion source which appear to come from a direction near the Sun and 
opposite the Sun, respectively. The next most active sources are the North and South toroidal which radiate from high above 
the ecliptic plane, about +60 degrees centered on the direction of motion about the Sun. The last two sources are the weakest 
in contribution but are the fastest particles, called the North and South Apex and are located +20 degrees also centered on the 
direction of motion about the Sun. Figure 4-11 depicts the radiant distribution at Earth of sporadic meteoroids for the larger 
threat size particles with masses greater than 100 micrograms. This figure is what the Earth sees on its trajectory around the 
Sun. Even though Mars is on a slightly more eccentric path the sources still appear to radiate from the same direction relative 
to the Sun and Mars’ motion. 
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Figure 4-11. Radiant distribution of sporadic meteoroids as seen by Earth in heliocentric ecliptic plane. Image from 
(Taylor and McBride, 1997). 


All environment files describe the cumulative cross-sectional flux of meteoroids down to a limiting mass of 10° grams. 
Particle density is assumed to be 1.0 g/cm”, the average density of cometary material, which comprises the majority of 
meteoric material in the inner solar system. 


The Meteoroid Environment Office specified the following relation for scaling the flux of 1 microgram meteoroids to larger 
meteoroid masses. 


Por > - (1) 


The variable F is the cumulative cross sectional area flux of meteoroids with mass m or larger. The parameter Fp 1s the 
cumulative cross sectional area flux of 1 microgram or larger meteoroids read from the environment files. Both F and Fo have 
units of number of meteoroids per square meter per year. The mass m has units of grams. 


The following subsections describe the method used to calculate the MEM meteoroid environment files for the mission 
phases listed in Table 4-6. 


4.4.3.1 From Earth to Mars Cruise Phase: 


The Interplanetary Meteoroid Engineering Model (IPMEM v 2.1) was used to calculate an average meteoroid environment 
file over the interplanetary cruise phase for Spirit, Mars Reconnaissance Orbiter (MRO) and a simple Hohman transfer. These 
three missions represent different trajectories to Mars with varying transfer times. State vectors were selected using the same 
time step between each of their Earth departures and Mars arrivals. These three trajectories were analyzed separately to 
generate an average meteoroid environment over the cruise portion of each mission, then the files were averaged to produce 
the environment file. 


Figure 4-12 is a plot of the cumulative cross sectional area flux of meteoroid one microgram or larger factored into 5 km/s 
closing speed bins for the outbound interplanetary cruise portion of the mission. 


23 


JSC- 66287 


Average Speed Distribution - Outbound 
Earth to Mars Interplanetary Cruise Phase 


8.0E-01 

= 

> 7.0E-01 

E 6.0€-01 

5 5.0€-01 

Lin 

@ 4.0E-01 

S 

= 3.0E-01 

a 

Y 2.0E-01 

» 

S 10€-01 

Ls 
0.0£+00 


0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 


Speedkm/s 


Figure 4-12. Cumulative flux of sporadic meteoroids distributed by 5 km/s speed bins for masses of one microgram or 
larger over the outbound interplanetary cruise portion of the mission; generated using average environments from 
the interplanetary cruises of Spirit, MRO and a simple Hohman transfer. 


4,4,3.2 Mars Arrival and Orbit 


There are 11 meteoroid environment files that describe the meteoroid environment during the aerobraking phase and final 
rendezvous phase before trans-Earth insertion. These files were generated using a newly developed meteoroid environment 
model, MarsMEM. This version of MEM includes Mars as the central body so that the correct planetary gravitational 
focusing and shielding will be used. The names of the environment files generated with MarsMEM are listed by mission 
phase in Table 4-6. The aerobraking phase of the trajectory was assumed to occur over 150 days. Ten orbits were constructed 
to simulate the falling apoapse. The Phase 12 file represents the final orbit where the vehicle is presumed to spend the most 
time before trans-Earth insertion. Each file is taken as representative of an “average” environment for that phase. 


Figure 4-13 is a plot of the cumulative cross sectional area flux of meteoroids | microgram or larger factored into 5 km/s 
closing speed bins for the Mars orbit phase of the mission. 
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Figure 4-13. Cumulative flux of sporadic meteoroids distributed by 5 km/s speed bin for masses of one microgram or 
larger onto a spacecraft in Mars orbit. 


4.4.3.3, From Mars to Earth Cruise Phase: 


For the interplanetary return phase of the mission, there were no specific trajectories described, and no recent NASA Mars 
missions to model this. A simple Hohman transfer was again used to describe the state vectors over the return trip. States 
were sampled in | day time increments over 8.5 months. 


Figure 4-14 is a plot of the cumulative cross sectional area flux of meteoroids | microgram or larger factored into 5 km/s 
closing speed bins for the Mars to Earth transit phase of the mission. 
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Figure 4-14, Average Speed Distribution over the inbound interplanetary cruise portion of the mission; generated 
using an average environment from the interplanetary cruises of a simple Hohman transfer orbit. 
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4,4,3.4 Meteoroid Environment Uncertainty 


Currently there is insufficient data to construct an “uncertainty” for the meteoroid environment. The one sigma uncertainty 
files that were provided by the Meteoroid Environment Office describe the environmental differences between the orbits or 
flight paths and do not reflect actual uncertainty in the environment. 


Earth to Mars transit phase — The Interplanetary Meteoroid Engineering Model (IPMEM v 2.1) was used to calculate an 
average meteoroid environment file over the interplanetary cruise phase for Spirit, Mars Reconnaissance Orbiter (MRO) and 
a simple Hohman transfer. These three missions represent different trajectories to Mars with varying transfer times. State 
vectors were selected using the same time step between each of their Earth departures and Mars arrivals. These three 
trajectories were analyzed separately and then averaged to produce the environment file. A 1 sigma variation was determined 
between the three different courses. The | sigma file can be added to or subtracted from the average interplanetary meteoroid 
file above to represent the variation in the meteoroid flux due to the unknown flight path 


Aerobraking and Mars orbit phases — The Aerobraking and the Mars orbit one sigma uncertainty files were generated 
using the customary MEM procedure. MEM is a Monte Carlo code that, in this case, generated a flux estimate for 60 second 
segments of the orbit. The various segments were averaged together to obtain an average flux for the trajectory. The standard 
deviation was calculated from the same data set. All orbits have the same inclination and approximately the same node since 
no other information was given. The MEO also said that the one sigma files were meant to bound the variation in flux due to 
unknown orbit plane orientation. 


Mars to Earth transit phase — The same uncertainty that was calculated from the outbound interplanetary cruise portion 
was used to develop the inbound sigma file. The one sigma standard deviation file for the Mars to Earth transit phase is 
intended to represent the variation in the meteoroid flux due to the unknown flight path. 


4.5 Ballistic Limit Equations 


4.5.1 SLA561V ablator BLE 


The SLA-561V impact tests were performed on two specimens provided by the Stardust program. The ablator was 0.75 
inches thick and was bonded to a 1.25 inch thick aluminum honeycomb. Specimen fabrication started with applying an 
adhesive film to the aluminum honeycomb face sheet. HRP F35-2.5 Flexcore was taped to the panels with J414 polyethylene 
tape to prevent it from moving. The assembly was then cured in an oven under an 8-inch Hg vacuum. The technicians then 
mixed a “wetcoat” of silicone resin, catalyst and solvent and the Super-Light Ablator SLA-561V. (The SLA material is 
composed of a silicone resin, catalyst, glass and cork fillers and has a bulk density of 0.257 g/cm” when cured). The test 
panels were then sprayed with DC-1200 silane primer. The wet coat was then sprayed onto the test panels. The mixed Super- 
Light Ablator (SLA-561V) was hand applied to fill the Flexcore cells. The test articles were then vacuum bagged, and heat 
cured. The tops of the core cells were machined flat. 


The test results are listed in Table 4-7. The first column lists the projectile and target materials. The second column lists the 
projectile diameter in cm. The third column lists the impact speed in km/s. The fourth column list the impact angle in degrees 
measured from the target normal. The fifth column lists the measured crater depth in cm. The sixth column lists the impact 
speed pi-group, the seventh column lists the crater depth pi-group. The last column lists the shot number. 


Table 4-7. SLA561V Hypervelocity Impact Test Results. 
d U B P d/"8 §°° (v cosB)”? = P/d Shot no. 
(cm) (km/s) (deg) (cm) 
Nylon->SLA 0.318 6.00 70 1.905 3.536322 6.0 A-2552 
2017-T4->SLA 0.060 6.83 0 0.670 10.12642 11.2 HITFO7001 


2017-T4->SLA 0.100 6.81 60 0.784 6.550056 HITFO7002 
2017-T4->SLA 0.100 6.88 0.928 8.309009 9.3 HITFO7003 
2017-T4->SLA 0.110 6.93 0 1.694 10.57520 15.4 HITFO7004 





Test A-2552 was performed in 1995 and the results are documented in (Gabriel and Christiansen, 1996). Four tests were 
performed on the second test article in 2007. The results of the second series of tests were document in (Lyons, 2007). Figure 
4-15 is a photograph of the second test article after the Orion program had performed 4 hypervelocity impact tests. 
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Figure 4-15. Second SLA561V Test Article at the Completion of Testing. 


The 1995 test point was fitted with the equation 
P 
a a d1/18 §°5 (V cosp)?/9 


reported in (Lyons, 2007), where P is the crater depth in cm, d is the projectile diameter in cm, 6 is the projectile diameter in 
g/cm’, V is the impact speed in km/s and is the impact angle measured from the target normal. 


This equation was coded up in BUMPER during November 2006 and was used for the aft TPS penetration analysis. The BLE 
and the test data are cross-plotted in Figure 4-16. 


27 


JSC- 66287 


PO 
qn 


O 2017-14 projectile 
© nylon projectile 
—fit 


ho 
CO 


— 
qn 


10 


penetration depth/projectile diameter, P/d 





0 2 4 6 8 
d1/18 §°.5 (V cosp)28 


Figure 4-16. Scaled SLA561V Crater Depth as a Function of Scaled Impact Speed. 


The failure criterion for the aft TPS was assumed to be a crater with depth equal to the 1.0 cm thickness of the aft TPS 
SLAS561V. The diameter of the aluminum orbital debris particle that will penetrate to 1.0 cm is plotted in Figure 4-17. 
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Figure 4-17. Orbital Debris Ballistic Limit Curve for 1.0 cm Deep Craters in Semi-Infinite SLA561V Targets. 
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Figure 4-18 Meteoroid Ballistic Limit Curve for 1.0 cm Deep Craters in Semi-Infinite SLA561V Targets. 


4.5.2 Carbon Phenolic BLE 


The carbon phenolic BLE was developed for this report from a Naval Research Laboratory report (Halperson, 1968), 
available from the Defense Technical Information Center, and from two tests performed by NASA for the Orion program 
(Hedman, 2007). 


Material characterization data were not readily apparent in the NRL reports. What was found is reproduced in Table 4-8. 


Table 4-8. Target and projectile material properties. 
Density UTS Sound speed 
(g/cm*) (KSI) (km/s) 

Carbon Phenolic, NRL 1.49 15 3.7 





Carbon Phenolic, Orion 1.48 3 
glass 2.50 


The glass projectile mass density was calculated from the nominal 1/8 inch diameter of the glass sphere projectiles and the 
reported projectile masses. (Variation in the nominal diameter of the spheres resulted in calculated densities ranging from 2.2 
to 2.7 g/cm’, with the average around 2.5 g/cm®. A 2.5 g/cm’ density indicates the projectiles were made of soda lime glass.) 
Reports prior to the one referenced mentioned testing carbon tape phenolic and chopped carbon tape phenolic. No mention 
was made in the referenced report which material was tested. However, a prior report mentioned that chopping the tape had 
more of an effect on the front surface spall than on the through hole diameter(Anon., 1963). 


A search of MATWEB was performed as a check on the Carbon Phenolic material properties listed in Table 4-8. The search! 
obtained a range of mass densities from 1.31 to 1.44 g/cm’ and a range of ultimate tensile strengths from 7.3 to 19 ksi, hence 
the values listed in Table 4-8 are within the range of possibilities. 


' http://www.matweb.com/search/DataSheet.aspx ?MatG UID=bdd98b907c6f4695b657 8b | 2efbd 1321 &ckck=1 
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The Carbon Phenolic test data extracted from the NRL report and the JSC test data are reproduced in Table 4-9. The first 
column lists the projectile and target materials. The second column lists the projectile diameter in cm, the third column lists 
the impact speed in km/s, the fourth column lists the crater depth in cm, the fifth column lists the impact angle in degrees 
measured from the target normal. The sixth and seventh columns list the pi-groups used to scale the impact speed and the 
crater depth respectively. The eighth column lists the page or figure number were the data were taken from in the NRL report. 


Table 4-9. 1/8 inch Glass and 3/16 inch Aluminum Spheres Impacting Carbon Phenolic. 


Projectile- D T y B a 
ee Ko1 (cm) (km/s) (cm) —_ (deg) V/V(¥/p) (p/8) P/d YelV Ig a= 


Glass->CP 0.318 085 0.143591 0 4.706248 0.451545 scaled from Figure 6 
Glass->CP 0.318 4.63 0.861549 0 25.57744 2.709273 scaled from Figure 6 


Glass->CP 0.318 6.18 0.963574 0 34.17146 3.030108 scaled from Figure 6 
Glass->CP 0.318 666 0.967353 0 36.83151 3.041991 scaled from Figure 6 


Glass->CP 0.318 8.03 1.073157 0 44.40243 3.374708 scaled from Figure 6 


Glass->CP 0.318 818 1.031591 0 45.22091 3.243998 scaled from Figure 6 
Glass->CP 0.318 5.76 0.842655 45 22.50814 2.649859 _ scaled from Figure 6 
2017->CP 0.477 6.81 1.09400 45 27.42584 2.293501 HITF10172 


2017->CP 0.421 6.77 1.07000 30 40.89712 2.541568 HITF10173 





The extra parameters in the column 6 pi-group are the target mass density p and the projectile mass density 6. 


The last two lines of Table 4-9 are the JSC data. These data were collected in 2010 and use an Orion Carbon Phenolic 
candidate material. The Carbon Phenolic was described in the test plan as: 


A structural carbon/phenolic composite pad provided by NASA Langley Research Center. The targets are 
machined to a 1” thick disk that is 9.5” in diameter. The carbon/phenolic pad has a density of approximately 92.3 
lb/ft? (1.48 g/cm”). The ply normal vector has an approximate 20° orientation with respect to the disk normal 
vector. The through-thickness Young’s modulus, shear modulus and tensile strength are approximately 2.6 Msi (18 
GPa), 0.9 Msi (6 GPa) and 3 ksi (20 MPa), respectively. 


The ultimate tensile strength (UTS) of the Langley Carbon Phenolic given above is 1/5 the UTS of the NRL material. 
The JSC test article post test is shown in Figure 4-19.The impacts at the 3 O’clock and the 9 O’clock positions were aimed 


near the edge to cause edge breakout. These tests are not reported here. The tests at the 12 O’clock and the 6 O’clock 
positions are HITF1072 and HITF1073, respectively. 
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Figure 4-19. Post Test Views of the JSC Test Article. 


Rear free surface spall in metals usually means that the crater depth in the finite thickness target is larger than the crater depth 
in a semi-infinite target for similar impact conditions. 


A close-up of the HITF1072 impact side crater is shown in Figure 4-20. 





Figure 4-20. HITF 1072 Impact Side Crater. 


The Carbon Phenolic data were scaled using pi-groups developed for metals cratering. The results are plotted in Figure 4-21. 
The metals scaling relation is plotted as the dashed curve. This curve applies to the crater depth in metals in the large impact 
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speed limit. There are reasons for thinking that Carbon Phenolic will follow the metals scaling relation at large impact 
speeds, however, the shock impedance mismatch between the projectile and target materials could preclude this from 
happening in the testable range of impact speeds. For aluminum projectiles striking aluminum targets the large speed limit 
the dashed curve applies to is 3 km/s and larger. However, if there is a large shock impedance mismatch between the 
projectile and the target, then the large speed limit could be at even larger speeds. For example, for nylon striking aluminum 
targets the large speed limit is 10 km/s and larger and for steel projectiles striking aluminum the large speed limit is 80 km/s 
and larger. The shock impedance of aluminum and glass is larger than the impedance of Carbon Phenolic so the authors 
expect the Carbon Phenolic hypervelocity impact data to lie above the metals crater depth scaling relation, which the NRL 
data do. However, the JSC data do not. This might indicate that for impact speeds less than the large speed limit the crater 
depth does not scale with the normal component of the impact velocity, like it does in the large speed limit. It might also 
indicate that there are material differences that are not captured by scaling with the UTS and the mass density. Further testing 
is required to establish a reliable method for extrapolating the crater depth data collected at testable impact speeds to impact 
speeds typical of meteoroids. 


For this preliminary assessment of the risk of forward TPS perforation the authors took the expedient of extrapolating the test 
data to large impact speeds using the metals crater depth scaling relation renormalized to lie above all the Carbon Phenolic 
test data. The solid curve in Figure 4-21 was the result. 
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Figure 4-21. Scaled Carbon Phenolic Crater Depth as a Function of Scaled Impact Speed. 
The solid curve in Figure 4-21 can be simplified to the following equation for Carbon Phenolic, 


P 
A 0.60885 d¥/18 §1/2(V cos B)?/3 


where P is the crater depth in cm, d is the projectile diameter in cm, 5 is the projectile diameter in g/cm’, V is the impact 


speed in km/s and f is the impact angle measured from the target normal. The above equation applies to the high strength 
material properties listed in Table 4-8 for the NRL carbon phenolic. 


The failure criterion for the forward TPS was assumed to be a crater with depth equal to the 1.2 cm thickness of the forward 
TPS Carbon Phenolic. The diameter of the aluminum orbital debris particle that will penetrate to 1.2 cm is plotted in Figure 
4-22. The SLAS561V BLE, plotted in Figure 4-17, and the Carbon Phenolic BLE, plotted in Figure 4-22, are plotted on the 
same scale. About a 6.5 times larger orbital debris particle 1s required to fail the forward TPS compared to the aft TPS, at 7 
km/s impact speed. The small size particles that lead to aft TPS failure are the reason why the aft TPS failure dominates the 
MMOD reliability of the MSREEV. 
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Figure 4-22. Orbital Debris Ballistic Limit Curve for 1.2 cm Deep Craters in Semi-Infinite Carbon Phenolic Targets. 
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Figure 4-23 Meteoroid Ballistic Limit Curve for 1.2 cm Deep Craters in Semi-Infinite Carbon Phenolic Targets. 
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4.5.3. Multi-Shock Shield BLE 


The MMOD “garage” shield shown in Figure 4-5 was analyzed with an advanced Nextel Multi-Shock shield that was 
developed at the NASA Johnson Space Center Hypervelocity Impact Technology Facility (HIT-F) (Christiansen, 1993). Of 
the three available configurations, only #2 was available in the Bumper code at the time of the analysis: 


1. Four equally spaced ceramic fabric bumpers with a flexible rear wall 
2. Four equally spaced ceramic bumpers with an aluminum rear wall 
3. Two equally spaced ceramic bumpers with a two-sheet aluminum Whipple shield 


The MMOD Shield Ballistic Limit Analysis program (Ryan and Christiansen, 2010) was used to size a shield that would 
protect against a 10 mm diameter meteoroid particle with a density of 1 g/cc and a velocity of 25 km/s impacting at 0 
degrees. The shielding design assumed a total shield stand-off distance of 30 cm, with three Nextel layers (total areal density 
= 0.19 g/cm’) equally spaced at 7.5 cm. The rear wall of the shield was a 0.18 cm 2219-T87 aluminum sheet (areal density = 
0.51 g/cm’) at a distance of 7.5 cm from the closest Nextel layer. The total areal density of this design was 0.698 g/cm* 
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Figure 4-24. 0° and 45° Ballistic Limit Curves for Nextel Multi Shock “Garage” Shield Concept sized to protect 
against a 1 cm meteoroid projectile at 25 km/s. Shield failure (detached spall) is predicted for particle diameters above 
the curves. 
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Three different studies were performed. First was a calculation of the risk of MMOD impact as a function of impact speed 
and impact angle. Second was a calculation of the probability of TPS failure. Third, and last, was the evaluation of an 


MMOD shield concept for protecting the EEV from meteoroid impact. 


5.1 MSR EEV Impact Analysis Results 


The impact analysis task produced tables of fluxes for discrete ranges of impact speed and impact angle. The impact speed is 
indexed by the row number and the impact angle is indexed by the column number. Tables were calculated for the forward 
and aft TPS, all 24 attitudes, mission phases 2 through 14, and four meteoroid diameters, for a total of 2,496 tables. Also, 
each mission phase was averaged over all 24 attitudes to produce an additional 104 tables. The four meteoroid diameters used 
to construct the tables are 1 mm, 100 microns, 50 microns, and 10 microns. 


Table 5-1 is the table of fluxes of meteoroids | mm or larger impacting on the forward TPS in Mars orbit. The first column of 
numbers is the row number. For the meteoroid tables, the first row corresponds to impact speeds from 0 to 5 km/s, the second 
to impact speeds from 5 to 10 km/s, the third to 10 to 15 km/s, ... to 95 to 100 km/s. The second column of numbers are the 
numbers of impacts onto the forward TPS during | year which strike the forward TPS at angles between 0 to 10 degrees 
measured with respect to the TPS normal. The third column is for impacts between 10 and 20 degrees, the fourth for impacts 
between 20 to 30 degrees, etc. 


Table 5-1 Table of Fluxes of 1 mm or Larger Meteoroids on the Forward TPS During the Mars Orbit Phase. See the 
text for a description of the table format. 
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One has to multiply Table 5-1 by the duration of exposure in years listed in Table 4-3 for the Mars orbit mission phase to 
obtain a number of impacts. Table 4-3 lists 3 years of exposure for the Mars Orbit phase, so one must multiply the fluxes in 
Table 4-3 by 3 to obtain the numbers of impact. If the Mars orbit phase is lengthened or shortened, then one would multiply 
by the new exposure duration. 


The first ten rows of Table 5-1 are plotted as a 3D column chart in Figure 5-1. These plots are called V-B plots at JSC/KX 
because impact speed is usually denoted as V and impact angle as B. The plot shows that randomizing the EEV attitude has 


smeared out the number of impacts over a broad range of impact angles. 
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Figure 5-1 v-B Plot for the Forward TPS in Mars Orbit. 


The tables for the flux of orbital debris particles were treated differently. First only one attitude was used, the 20 degrees 
entry angle shown in Figure 4-7. Second only one mission phase was analyzed, mission phase 15 Earth entry. Mission phase 
1, launch to LEO, was not analyzed. This gives a total of 4 tables for the orbital debris analysis of the forward TPS, one for 
the 1 mm OD particles, one for the 100 micron diameter particles, one for the 50 micron diameter particles and one for the 10 
micron diameter particles, and 4 tables for the aft TPS, for a total of 8 tables. 


The orbital debris tables also differ from the meteoroid tables by the speed increment used for each row. The orbital debris 
tables have 20 rows, just like the meteoroid tables, and the row number is listed in the first column, just like the meteoroid 
tables, but row | corresponds to 0 to 1 km/s, row 2 to | to 2 km/s, ..., row 20 to 19 to 20 km/s 


The average attitude V-B tables for each mission phase were multiplied the duration of exposure listed in Table 4-3, summed 
Over impact speed and angle, summed over both sides of the EEV, and then summed over the 15 mission phases to obtain a 
total number of impacts for each of the 4 particle sizes analyzed. The results are shown in Table 5-2 broken out by orbital 
debris particle impacts and by meteoroid impacts. 
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Attitude 


a" 


WO CON DW BPW N 


6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 
6.78E-06 


6.78E-06 
6.78E-06 


Table 5-2. MSR EEV Analysis - MMOD Impact Study Results Summary. 


7.54E-04 
7.13E-04 
4.97E-04 
5.62E-04 
4.90E-04 
6.70E-04 
6.93E-04 
5.57E-04 
6.34E-04 
6.87E-04 
7.19E-04 
7.29E-04 
5.44E-04 
5.75E-04 
5.11E-04 
5.67E-04 
6.39E-04 
4.77E-04 
5.76E-04 
7.94E-04 
7.54E-04 
4.87E-04 


5.20E-04 
7.28E-04 


Number of Impacts - per Limiting Particle Size 


0.001 cm 


7.60E-04 
7.20E-04 
5.03E-04 
5.69E-04 
4.96E-04 
6.76E-04 
6.99E-04 
5.64E-04 
6.40E-04 
6.94E-04 
7.26E-04 
7.35E-04 
5.51E-04 
5.81E-04 
5.17E-04 
5.74E-04 
6.46E-04 
4.84E-04 
5.83E-04 
8.01E-04 
7.60E-04 
4.93E-04 


5.27E-04 
7.35E-04 


OD 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 
9.00E-05 


9.00E-05 
9.00E-05 


MM 
1.22E-02 
1.16E-02 
8.06E-03 
9.12E-03 
7.94E-03 
1.09E-02 
1.12E-02 
9.04E-03 
1.03E-02 
1.11E-02 
1.17E-02 
1.18E-02 
8.82E-03 
9.32E-03 
8.28E-03 
9.20E-03 
1.04E-02 
7.74E-03 
9.35E-03 
1.29E-02 
1.22E-02 
7.89E-03 


8.44E-03 
1.18E-02 


Total 
1.23E-02 
1.17E-02 
8.15E-03 
9.21E-03 
8.03E-03 
1.10E-02 
1.13E-02 
9.13E-03 
1.04E-02 
1.12E-02 
1.18E-02 
1.19E-02 
8.91E-03 
9.41E-03 
8.37E-03 
9.29E-03 
1.05E-02 
7.83E-03 
9.44E-03 
1.30E-02 
1.23E-02 
7.98E-03 


8.53E-03 
1.19E-02 


OD 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 
5.40E-03 


5.40E-03 
5.40E-03 


MM 
7.89E+00 
7.47E+00 
5.20E+00 
5.89E+00 
5.13E+00 
7.01E+00 
7.25E+00 
5.83E+00 
6.63E+00 
7.19E+00 
7.53E+00 
7.63E+00 
5.70E+00 
6.02E+00 
5.35E+00 
5.94E+00 
6.70E+00 
5.00E+00 
6.04E+00 
8.31E+00 
7.89E+00 
5.10E+00 


5.45E+00 
7.62E+00 


Total 
7.90E+00 
7.47E+00 
5.21E+00 
5.89E+00 
5.13E+00 
7.02E+00 
7.26E+00 
5.84E+00 
6.64E+00 
7.20E+00 
7.53E+00 
7.64E+00 
5.70E+00 
6.02E+00 
5.35E+00 
5.94E+00 
6.70E+00 
5.00E+00 
6.04E+00 
8.32E+00 
7.90E+00 
5.10E+00 


5.45E+00 
7.63E+00 


OD 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 
1.24E-01 


1.24E-01 
1.24E-01 


MM 
8.26E+04 
7.82E+04 
5.45E+04 
6.16E+04 
5.37E+04 
7.34E+04 
7.60E+04 
6.11E+04 
6.95E+04 
7.53E+04 
7.88E+04 
7.99E+04 
5.96E+04 
6.30E+04 
5.60E+04 
6.22E+04 
7.01E+04 
5.23E+04 
6.32E+04 
8.70E+04 
8.26E+04 
5.34E+04 


5.70E+04 
7.98E+04 
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Total 
8.26E+04 
7.82E+04 
5.45E+04 
6.17E+04 
5.37E+04 
7.34E+04 
7.60E+04 
6.11E+04 
6.95E+04 
7.53E+04 
7.88E+04 
7.99E+04 
5.96E+04 
6.30E+04 
5.60E+04 
6.22E+04 
7.01E+04 
5.23E+04 
6.32E+04 
8. 70E+04 
8.26E+04 
5.34E+04 


5.70E+04 
7.98E+04 


asco) [amin (500gm) | mm at0ym) [tm 0pm) 
| max | 6.78E-06 _7.94E-04 8.01E-04 | 9.00E-05 1.29E-02 __1.30E-02 
9.00E-05 1016-02 __1.016-02 
| min | 6.786-06 __ 4.776-04_4.846-04] 9.00E-05 _7.74E-03__7.83E-03 


The second from the bottom row lists the number of impacts averaged over all 24 attitudes and all 15 mission phases, 1.e., the 
average of the column of 24 numbers above it. The third row from the bottom lists the maximum of the column of 24 
numbers and is representative of the mean number of impacts when in a “worst” case attitude. The last row lists the minimum 
of the column of 24 numbers and is representative of the mean number of impacts when in a “best” case attitude. 





Note that the number of orbital debris impacts relative to the number of meteoroid impacts increases with the size of the 
particle. The 10 micron orbital debris particles impacts are 0.000119% of the 10 micron meteoroid impacts, 0.0684% of the 
50 micron size particle impacts, 0.57% of the 100 micron size particle impacts and 0.63% of the 1 mm size particle impacts. 
However, in every case it is a significantly smaller number than the number of meteoroid impacts. 


5.1.1. 1.0 mm Impact Results 


The results for the MSR EEV impact study for a limiting particle size of 1.0 mm are shown in Table 5-3. The results are 
broken out for each mission phase as well as for the two MSR EEV TPS components (aft and forward). 


Table 5-3. MSR EEV Analysis - Impact Study Results for 1.0 mm Limiting Particle Size. 


MSR EEV Analysis — 1.0 mm Impact Results 





Number of Impacts - per Limiting Particle Size 
Missions Phase | Forward | Aft | Total 


Launch-LEO (not analyzed) 


Z E arthtoMarsTransit_—~S*~<~*~dSCS*é~‘iC CE 3986-05 ~~ T69E-04 
a Mars Aeropra ing merlot QZE05 | 9.O7E-06 | 3.82E-05 
2 Mars Orbit 938 E04 |S 99E-OS | 2.53E-04 
i MarstoEarth Insert = | A OZE-O7 | 125E-07 | 5 27E-O7 
Mars to Earth Transit 1.22E-04 3.74E-05 1.59E-04 
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MSR EEV Analysis — 1.0 mm Impact Results 
Number of Impacts - per Limiting Particle Size 
Missions Phase | Forward | Aft =| Total 


EEV Entry (OD only) 3.13E-06 3.66E-06 6.78E-06 
4.77E-04 1.50E-04 6.26E-04 





5.1.2 0.5 mm Impact Results 


The results for the MSR EEV impact study for a limiting particle size of 0.5 mm are shown in Table 5-3. The results are 
broken out for each mission phase as well as for the two MSR EEV TPS components (aft and forward). 


Table 5-4. MSR EEV Analysis - Impact Study Results for 0.5 mm Limiting Particle Size. 
MSR EEV Analysis — 0.5 mm Impact Results 


Number of Impacts - per Limiting Particle Size 
Missions Phase | Forward | Aft | Total 


Launch-LEO (not analyzed) 





Earth to Mars Transit 
Mars Aerobraking 
Mars Orbit 

Mars to Earth Insert 
Mars to Earth Transit 


2, 10E- 03 
4.73E-04 
3.13E-03 
6.53E-06 
1.97E-03 


6. ASE. 04 
1.47E-04 
9.72E-04 
2.03E-06 
6.07E-04 
4.55E-05 


2. 74E- 03 
6.20E-04 
4.10E-03 
8.56E-06 
2.58E-03 
9.00E-05 


EEV Entry (OD only) 4.46E-05 





5.1.3. 0.1 mm Impact Results 


The results for the MSR EEV impact study for a limiting particle size of 0.1 mm are shown in Table 5-5. The results are 
broken out for each mission phase as well as for the two MSR EEV TPS components (aft and forward). 


Table 5-5. MSR EEV Analysis - Impact Study Results for 0.1 mm Limiting Particle Size. 
MSR EEV Analysis — 0.1 mm Impact Results 


Number of Impacts - per Limiting Particle Size 
Missions Phase | Forward | Aft =| Total 


Launch-LEO (not analyzed) 

Earth to Mars Transit 1. 35E+00 4. 16E. 01 i 77E+00 
Mars Aerobraking 3.05E-01 9.49E-02 4.00E-01 
Mars Orbit 2.02E+00 6.27E-01 2.65E+00 
Mars to Earth Insert 4.21E-03 1.31E-03 5.52E-03 
Mars to Earth Transit 1.27E+00 3.92E-01 1.67E+00 
EEV Entry (OD only) 3.09E-03 2.31E-03 5.40E-03 


4.95E+00 1.53E+00 6.50E+00 


5.1.4 0.01 mm Impact Results 


The results for the MSR EEV impact study for a limiting particle size of 0.01 mm are shown in Table 5-6. The results are 
broken out for each mission phase as well as for the two MSR EEV TPS components (aft and forward). 


Table 5-6. MSR EEV Analysis - Impact Study Results for 0.01 mm Limiting Particle Size. 
MSR EEV Analysis — 0.01 mm Impact Results 


Number of Impacts - per Limiting Particle Size 
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Launch-LEO (not analyzed) se sue Loe 
Earth to Mars Transit 1.42E+04 4.36E+03 1.85E+04 
Mars Aerobraking 3.20E+03 9.94E+02 4.19E+03 


Mars Orbit 2.12E+04 6.57E+03 2.77E+04 


Mars to Earth Insert 4.41E+01 1.37E+01 
Mars to Earth Transit 1.33E+04 4.10E+03 1.74E+04 
FEV Entry (OD only) 7.25E-02 5.17E-02 





5.1.5 Metrics 


The probability of meteoroid impact for a given diameter meteoroid or larger is plotted in Figure 5-2. The diameter meteoroid 
that impacts with a probability equal to the requirement, when the EEV is the lower risk attitude is 5.8 mm, and the diameter 
meteoroid that strikes with a probability equal to the higher risk attitude is 6.6 mm. 


1.E+00 


—low risk (attitude 18) 
— high risk (attitude 20) 


— -requirement 


1.E-01 


1.E-02 


1.E-03 


1.E-04 


1.E-05 


0.661 cm 


probability of impact 


15-06 |} 0.583 cm 


1.E-07 





1.E-08 
0.001 0.01 0.1 1 10 
micrometeoroid diameter (cm) 


Figure 5-2 Diameter meteoroid that has a probability of 4x10" of impacting the EEV during 4.38 years. 


The simple meteoroid power-law flux relation Eq. (1) means that if the number of impacts N,e¢ at reference diameter d,.,, that 
the number of impacts by meteoroids with diameter d or larger is 


d —4.02 


The fraction of meteoroid impacts that occur during mission phase are plotted in Figure 5-3. The part due to impacts on the 
forward TPS are plotted in red and the part that are due to impacts on the aft TPS are plotted in blue. The orbital debris 
cumulative flux relation is not a power-law, so number of orbital debris impacts cannot be cross plotted with the meteoroid 
results in Figure 5-3. 
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Figure 5-3 Fraction of meteoroid impacts at a given size as a function of mission phase and whether the impact was on 
the forward or aft TPS. 


5.2 MSR EEV Probability of TPS Failure by MMOD Impact 


Table 5-7. is a list of the results from a calculation of the probability of TPS failure by MMOD impact. The first column lists 
the mission phase, the second column lists the probability of an MMOD impact onto the forward TPS causing a TPS failure 
on Earth entry, the third column lists the probability of an MMOD impact onto the aft TPS causing a TPS failure on Earth 
entry and the fourth column is the probability of an impact on either the forward or aft TPS causing a TPS failure. All three 
columns are probabilities calculated by averaging over all 24 attitudes listed in Figure 4-6. 


Table 5-7. MSR EEV Probability of TPS Failure Due to MMOD Impact. 
MSR EEV TPS Failure by MMOD Impact 


Probability of TPS Failure 
Missions Phase | Forward | Aft | Total 


Launch-LEO (not analyzed) a 

Earth to Mars Transit i. 326- 06 5. S36. 04 5. SAE- 04 
Mars Aerobraking -O5E- | 8.90E-05 8.92E-05 
Mars Orbit | .33E- | 5.78E-04 5.79E-04 
Mars to Earth Insert .74E- | 1.19E-06 1.19E-06 
Mars to Earth Transit 5.01E-04 5.03E-04 
au Entry (OD only 1 15E-06 4.09E-06 5.14E-06 


|“worst’caseattitude | | | BOE-O3 
best” caseattitude (| | iSCSI 
|Requirement | | | OOE-O7 

















The important take-away from Table 5-7 is that the probability of failure of the TPS due to MMOD impact exceeds the 
4x107' reliability allocation. But it is also significant that the probability of failure due to orbital debris impact during Earth 
entry is larger than the allocation. This is problematic because the TPS is not shadowed/shielded by the orbiter following 


40 


JSC- 66287 


separation. Separation times for prior spacecraft sample return canisters were approximately 4 hours for Genesis and Stardust 
and approximately 3 hrs for Hayabusa. In the case of Genesis, separation occurred outside of geosynchronous orbit and hence 
outside of the Earth orbital debris belt. If the EEV trajectory is similar, then the EEV will traverse the entire Earth orbital 
debris belt without the protection of the orbiter. If separation could be deferred to a later time, then the orbiter could provide 
more protection. The Apollo CSM separation occurred about 20 minutes before Earth interface. If the Genesis/SRC 
separation was deferred to 20 minutes before Earth interface, then separation would have occurred around 7,000 km. 
According to Figure 4-10, this would reduce the total fluence of 1 mm particles by a 1/3, which would help. 


Another take-away from Table 5-7 is that the probability of failure due to perforation of the aft TPS is 330 times larger than 
the probability of failure of the forward TPS. This is a consequence of the forward and aft TPS ballistic limit equations. The 
size meteoroid that penetrates the aft TPS at 7 km/s and normal incidence is 0.847 mm and the size meteoroid that penetrates 
the forward TPS at 7 km/s is 5.56 mm. Substituting the 0.847 mm diameter into the meteoroid flux equation results in a flux 
1,930 times larger than the flux of 5.56 mm diameter meteoroids, so it’s not surprising that there is a much greater probability 
of aft TPS failure than forward TPS failure. 


The line labeled “worst” case attitude in Table 5-8 is the attitude from Table 4-4 that produced the largest probability of TPS 
failure. The line labeled Total is the average of probability of TPS failure over the 24 attitudes listed in Table 4-4. This is 
representative of the probability of TPS failure for a randomly oriented EEV. The line labeled “best case attitude from Table 
4-4 is the attitude that produced the smallest probability of TPS failure. 


A bar chart similar to Figure 5-3 can be constructed for the mean number of penetrations (which is approximately equal to 
the probability of failure). The result is shown in Figure 5-4 the bar chart illustrates that the aft TPS is failing the requirement 
by 3 orders of magnitude. 
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1.E-04 





0.E+00 


launch-LEO = Earthto Mars Mars Orbit Mars Mars to EEV Entry 
(not Mars Transit Aerobraking Departure Earth Transit (OD) 
analyzed) orbit 


Figure 5-4 Mean Number of TPS Penetrations as Function of Mission Phase. 


5.3 MMOD Shield Concept 


The results for the MSR EEV MMOD shield concept PNP study are shown in Table 5-8. The results are broken out for each 
mission phase as well as for the two MSR EEV TPS components (aft and forward) and the conceptual MMOD shield. 
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Table 5-8. MSR EEV Analysis - MMOD Shield Concept Probability of TPS Failure due to MMOD Impact. 


MSR EEV MMOD Shield Concept Results 


Number of Perforations 
Missions Phase | Forward | At. =~ | ~—s MMOD Shield 
Launch-LEO (not analyzed) 


Earth to Mars Transit 0.00E+00 0. O0E+00 . 8. O1E- 08 8. O1E- 08 
Mars Aerobraking | 0. Econ | 0.00E+00 1.43E-08 | 1.43E-08 
Mars Orbit 0.00E+00 9.44E-08 9.44E-08 


Mars to Earth Insert ; 0.00E+00 | 1.97E-10 . 1.97E-10 
Mars to Earth Transit 0. OOE+00 0.00E+00 7.50E-08 7.50E-08 


— Entry (OD ont) 1.15E-06 4.09E-06 0.00E+00 5.24E-06 


|“worst’caseattitude | | | EOE 
Total | SEG | 4.09E-06 | 2.64E-07 |S SOE-06 
"best caseattiude | | SSSC~iSC“‘“‘SSCNWC#”«SOC“*@*‘COSE-OG 
|Requirement [POET 





The analysis was performed assuming that the probability of TPS failure given an MMOD shield perforation was one. This is 
a conservative assumption. The assumption is necessary because testing is required to determine the extent of TPS damage 
given an MMOD shield perforation. 


Table 5-8 shows that the MMOD shield has reduced the risk of TPS damage due to meteoroid impact to 2.64x10’, which is 
below the 4x10’ requirement. However, the possibility of orbital debris impact occurs after the EEV is released from the 
orbiter and the EEV is no longer protected by the MMOD shield. Thus the risk of TPS failure due to orbital debris impact is 
unchanged from the unshielded cased, section 5.2. This short 2 hour exposure to the orbital debris environment is sufficient 
to cause the EEV to fail its probability of TPS damage due to MMOD impact by a wide margin. 


Contours of numbers of TPS failures per square meter per exposure duration times 100 are plotted in Figure 5-5 for the 
orbital debris impacts during Earth entry phase. The EEV velocity vector is in the XZ plane of the figure and is pointing out 
the apex of the forward TPS cone. The contours of number of TPS failures are concentrated on the starboard side of the EEV 
due to the assumed argument of perigee of the orbit. The “hot” spot on the aft TPS will appear in a different location on the 
energy-absorbing sphere with different arguments of perigee. 


The contouring levels preclude seeing the impact contours on the forward TPS because the number of impacts on the forward 
TPS are a %4 of the number of impacts on the aft TPS. 
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Looking +X — , Looking -X 





Figure 5-5. MSR EEV Analysis —- MSR EEV Risk Contour Plot for Earth Reentry Phase 
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6. Conclusions 


The MSR EEV MMOD analysis has shown that the currently designed MSR EEV thermal protection system does not meet 
the 4x10’ maximum allowable risk requirement, given the assumptions made in this report. The aft TPS is made of 
SLA561V and fails the requirement by a factor of 4,300. The small density and strength of SLA561V makes an MMOD 
shield a necessity to meet the stringent 4x10” probability of TPS failure requirement. The forward TPS is made of carbon 
phenolic which is robust to MMOD impact, however, it too fails the requirement, but by the smaller factor of 13. 


The assumption that perforation of the TPS down to the base structure leads to loss of specimen containment needs to be 
confirmed by TPS designers. The Orion program assumes smaller depths of penetration into the TPS results in loss of 
vehicle, so complete penetration of the TPS is likely an optimistic assumption. The Orion program TPS designers provided a 
map of allowable TPS penetration depths based on thermal allowable of the tile bond lines. Something similar is needed for 
EEV. 


Surprisingly, the orbital debris impacts during entry was about 22% of the forward TPS risk of failure. This number may be 
an overestimate due to the conservative extrapolation of the orbital debris environment from 2,000 km to geosynchronous 
orbit made by Bjorkman. It may be possible to eliminate this conservatism by tasking the Orbital Debris Program Office with 
modifying their ORDEM 3.0 orbital debris environment to analyze the EEV Earth entry trajectory. ORDEM 3.0 does not 
require an extrapolation above 2,000 km altitude and greater confidence could be attached to the results. 


Finally, hypervelocity impact tests of the actual EEV TPS materials are recommended. If it is desired to optimize the weight 
of the MMOD shield for the EEV, then testing to determine the conditions leading to unacceptable TPS damage following 
perforation of the shield is needed. 


6.1 Forward Work Summary 
The following summarizes forward work needed to reduce uncertainty in the MMOD risk assessment. 


o MMOD requirements — need to confirm and/or update the 4x10°’ maximum risk value used in this assessment. 

o Spacecraft TPS — thickness and material type for the EEV TPS need to be confirmed and/or updated in the next 
MMOD assessment. 

o TPS failure criteria — need further efforts to define TPS failure criteria (damage limits for reentry). This will likely 
require tests (Ampact followed by arc-jet) and modeling of heating surrounding TPS cavities during reentry. 

o Assessment inputs — need to include analysis of meteoroid and orbital debris impacts during initial launch/LEO 
phase, and meteoroids during reentry phase. Even though the MMOD risk during these mission phases is believed to 
be small relative to other mission phases, it should be quantified. Need to update the analysis using ORDEM 3.0 
after it 1s released (no earlier than end of 2011). 

o Ballistic limit equations — additional tests and simulations should be considered to reduce uncertainty in the ballistic 
limit equations used in this analysis, particularly TPS response at impact velocities greater than 7 km/s. The ballistic 
limit equations for the aft TPS need to be reviewed and updated based on additional hypervelocity impact tests. 

o Optimize MMOD shielding/bioshield cover — evaluate additional shielding concepts to reduce mass and improve 
survivability. Evaluate risk as a function of when the MMOD shield is separated from the EEV at the end of 
mission. 

© Quantify uncertainty — assess uncertainty in risk assessment results using MEM 1-sigma files. 

o Other items described throughout the report and mentioned in conclusions that should be done to reduce uncertainty. 
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8. Appendix A - MSR EEV MMOD Impact Risk Analysis Statement of Work 


Statement of Work 
Analysis of Micrometeorite (MM) impact on the Mars Sample Return Earth Entry Vehicle (MSREEV) 


An agreement between NASA Ames Research Center 
March 7’ 2011 


JSC POC: Eric Christiansen/ ARC POC: Ethiraj Venkatatpathy 


Information to be provided by NASA Ames: 
1. Approximate MSR mission concept. 
Launch Date of the spacecraft delivering the MSREEV to Mars 
Arrival date of the MSREEV and its orbit and duration in orbit at Mars 
Departure date from Mars and the Date of Earth Entry 
Drawing of the MSREEV as conceived in the ~ 2000 time frame 
No information is currently available for the spacecraft that will deliver the MRSEEV to Mars nor that on 
which it will be returned 
to the Earth 


ono oR 


2. MSR EEV Cross Section with Materials lay-up 


Analysis to be provided by NASA JSC in a written report: 
1. Description of the MM distributions for the MSR missions and bounds on them 
MMOD in Low Earth Orbit (LEO) on way to Mars 
MM for transit to Mars 
MM during Mars Orbit 
MM during transit to Earth from Mars 
MMOD distribution in LEO upon return from Mars 


cnose 


2. Probabilities of hypervelocity impact as a function of velocity and particle sizes for the dwell times in each mission 
phase and the entire mission. 


Transit through LEO 

Transit from LEO to Mars Orbit 

While in Mars Orbit 

Transit from Mars to Earth 

Transit through LEO to Earth’s surface 

Integrated probability of MMOD impact as a function of velocity and particle size for the entire mission. 


monose 


3. Discussion of how the results of the analysis would change for missions of similar definition for later 
implementation (shifts on 26 month cycles) 


4. Discussion of MMOD impact testing on Carbon Phenolic and 2D carbon-Carbon from prior HVI testing. 


5. Discussion of directional aspects of the MM and MMOD fluxes and how the MSREEV could be best oriented with 
respect to the transit spacecraft(s) to avoid penetration of the outer mold line. 


Written report due 6 months after receipt of ISPT funds at JSC 


Interim presentation of results at 4 months from authority to proceed. 
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9. Appendix B - JPL Horizons Genesis Spacecraft Trajectory File Header 


KKKKKKKKKKKKRKKKKRKKKKRKKKKRKKKKKKKKKRKKKKRKKKKRKKKKRKKKKKRKKK KR KKK KR KKK KKKKKKKKKKRKKKKKKEKK 


Revised: Nov 19, 2004 Genesis (SRC) Spacecraft / (Earth) —47900 
http://genesismission.jpl.nasa.gov/ 


Genesis Sample Return Cannister (SRC): 
Launched: Aug 8, 2001 12:13:40 p.m. EDT (Delta 2 rocket) 


The SRC release was a one-second duration event beginning at 8-Sep-2004 
Tizs4:1il.eS Cl. End of Contact with the springs was to occur at 11*54:12.5. Cr. 
At this time, object -47900 (Genesis Sample Return Capsule) separated from 
-47 (Genesis bus). Half an hour later, the bus was diverted to an ultimately 
heliocentric trajectory, available separately using ID -47. 


SRC atmospheric entry was Sep 8, 2004 15:52:47.6 UTC. 

The SRC failed to deploy the drogue/parafoil at entry +127 seconds as planned, 
instead tumbling ballistically through the Earth's atmosphere to impact in 
sandy soil at geodetic coordinates 40.1278 deg (lat), 246.4919 deg (long). 


SPACECRAFT TRAJECTORY (from JPL navigation team): 


Solution name Span 
OD145 2004-Aug-09 to 2004-Sep-06 
OD154 2004-Sep-06 to 2004-Sep-08 


KKKKKKKKKKKKRKKKKRKKKKRKKKKRKKKKKKKKKRKKKKRKKKKRKKKKRKKKKKRKKK KR KKK KKK KKRKKKKKRKKKKKKKKKKEKK 


KKKKKKKKKKKKRKKKKRKKKKRKKKKRKKKKKKKKKRKKKKRKKKKRKKKKRKKKKKRKKK KR KKK KKK KKKKKKKRKKKKKKKKKKEKK 


Ephemeris / WWW_USER Mon Sep 12 12:09:52 2011 Pasadena, USA / Horizons 
KKEKKKKKKKKKRKKKKEKKKKKKKE KKK KKK KKK KKK KEK KKK KKK KKK KKK KEK KKK KKK KKK KK KEKE KE IK KKK KKK EK EK EK 
Target body name: Genesis (SRC) Spacecraft (-47900) {source: pfile_src_od154.nio} 
Center body name: Earth (399) {source: DE405} 

Center-site name: BODY CENTER 


KKKKKKKKKKKKRKKKKRKKKKRKKK KR KKK KKK KKKRKKKKRKRKKKRKKKKRKKKK KR KKK KR KKK KR KKK KKKKKKRKKKKRKKKKKKEKK 


Start time ¢ Ax Ds. 2004—=sep-07 13:52:00.0000 Cr 
Stop time ; AeDy ZO04—SeEp=-Us 15752:00.0000 -CTr 
SLEp=size >: 1 minutes 


KKKKKKKKKKKKRKKKKRKKKKR KKK KKK KKKKKKRKKKKRKKKKRKRKKKRKKKKKR KKK KKK KKRKRKKKKKKKKRKKKKKKKKKKEKK 


Center geodetic : 0.00000000,0.00000000,0.0000000 {E-lon (deg) ,Lat (deg) ,Alt (km) } 
Center cylindric: 0.00000000,0.00000000,0.0000000 {E-lon(deg),Dxy (km) ,Dz (km) } 


Center radii : S370. & Cote. & G350.0. kin {Equator, meridian, pole} 
OUCDUEL Units : KM=S 

Output format ; OS 

Reference frame : ICRF/J2000.0 

Output type >: GEOMETRIC cartesian states 


Coordinate systm: Earth Mean Equator and Equinox of Reference Epoch 


KKKKKKKKKKKKRKKKKRKKKKRKKKKRKKKKKRKKKKRKKKKRKKK KR KKK KRKKKKKR KKK KR KKK KKK KKRKKKK KR KKK KR KKKKKKEKK 


Coordinate system description: 


EFarth Mean Equator and Equinox of Reference Epoch 
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Reference epoch: J2000.0 


xy-plane: 
x=axis 


Z-ax1s 


plane of the Earth's mean equator at the reference epoch 

out along ascending node of instantaneous plane of the Earth's 
orbit and the Earth's mean equator at the reference epoch 
along the Barth mean north pole at the reference epoch 


Symbol meaning 


JDCT 
X 
oY 
Zi 
VX 
AG 
VZ 
LT 
RG 
RR 


Epoch Julian Date, Coordinate Time 
x-component of position vector (km) 
y-component of position vector (km) 
z-component of position vector (km) 
x-component of velocity vector (km/sec) 
y-component of velocity vector (km/sec) 

z-component of velocity vector (km/sec) 

One-way down-leg Newtonian light-time (sec) 

Range; distance from coordinate center (km) 
Range-rate; radial velocity wrt coord. center (km/sec) 


Geometric states/elements have no aberration corrections applied. 


Computations: by 
Solar System Dynamics Group, Horizons On-Line Ephemeris System 
4800 Oak Grove Drive, Jet Propulsion Laboratory 
Pasadena, CA 91109 USA 
Information: http://ssd.jpl.nasa.gov/ 


Connecr 


Author 


telnet://ssd.jpl.nasa.gov:6775 (via browser) 
telnet ssd.jpl.nasa.gov 6775 (via command-line) 
Jon.Giorgini@jpl.nasa.gov 


KEKKKKKKKRKKKKRKKKKRKKKKRKKKKKRKKKKRKKKKRKKKKRKRKKKKKKKKRKKKKRKKKKRKRKKKKKKKKRKKKKRKKKKKKKKKKK 
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10. Appendix C - SLA561V Ablator Ballistic Limit Equation 


Fortran code of the SLA561V ablator ballistic limit equation implemented in bumper: 


SLA-561V penetration equation (based on 1 data point, 

nylon projectile, 7Odeg impact angle, 6km/s impact into 
1.9cm thick SLA-561V bonded to honeycomb panel 

damage description indicated entry hole 2.4cm x 3cm in size, 
some SLA filled cells completely destroyed, others damaged, 
penetration to bondline assumed but may be less). 


Damage mode: Penetration into SLA-561V only (not into substrate) 


source: 
B.L.Christiansen, 13 Nov 2006. 


Coded By: 
Tom Prior, Hamilton Sundstrand/ESCG, 13-NOV-2006 
SUBROUTINE SLA561V(diam,velm,angr,projden, & 
& ABLATOR_THICK, ABLATOR_DEPTH) 


REAL (KIND=4) , INTENT (OUT) :: DIAM 

REAL (KIND=4) , INTENT(IN) :: VELM, ANGR, PROJUDEN 

REAL (KIND=4) ,INTENT(IN) :: ABLATOR_THICK, ABLATOR_DEPTH 
REAL (KIND=4) :: VNORM 


REAL*8 PDEPTH 
CALCULATE CRITICAL PENETRATION DEPTH 
PDEPTH = ABLATOR_THICK * (ABLATOR_DEPTH/100.0) 
Calculate the normal component of the velocity vector 
VNORM=VELM* COS (ANGR) 
CALCULATE CRITICAL PROJECTILE DIAMETER IN CM 


DIAM= (PDEPTH/3.7/SQRT (PROJDEN) /VNORM** (2./3.))**(18./19.) 
END SUBROUTINE SLA561V 
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11. Appendix D - Carbon Phenolic Ablator Ballistic Limit Equation 


> @brief Carbon Phenolic TPS BLE. 

! @details The equation was developed by Michael D. Bjorkman, JSC/ESCG. 

! The BLE derivation is documented in JSC-66287. It 1s based on 7 tests 

! performed at NRL with glass projectiles and 2 tests performed for JSC 

! using 2017-T4 projectiles. 

@param[out] diam projectile diameter in cm 

| @param[in] velm impact speed in km/s 

: @param[in] angr impact angle measured for the element normal in radians 
| (OG B17 2) 

| @param[in] projden projectile mass density in g/cm%*3 

! @param[in] ablator_thick carbon phenolic ablator thickness in cm 

| 
| 
| 


when penetrated causes failure. 
@section sec2 Change History: 
\li /AUG 8, 2011/ subroutine added by Michael D. Bjorkman 


SUBROUTINE CarbonPhenolic(diam,velm,angr,projden, & 
& ABLATOR_THICK, ABLATOR_DEPTH) 


REAL (KIND=4) , INTENT (OUT) :: DIAM 

REAL (KIND=4) , INTENT(IN) :: VELM, ANGR, PROJDEN 

REAL (KIND=4) ,INTENT(IN) :: ABLATOR_THICK, ABLATOR_DEPTH 
REAL (KIND=4) :: VNORM 

REAL (KIND=8) :: PDEPTH 


! CALCULATE CRITICAL PENETRATION DEPTH 
PDEPTH = ABLATOR_THICK * (ABLATOR_DEPTH/100.0_DP) 

!' Calculate the normal component of the velocity vector 
VNORM=VELM*COS (ANGR) 


! CALCULATE CRITICAL PROJECTILE DIAMETER IN CM 


@param[in] ablator_depth fraction, in percent, of the ablator thicknss that 


DIAM= (PDEPTH/0.60885_DP/SQRT (PROJDEN) /VNORM** (2.0_DP/3.0_DP))** (18.0_DP/19.0_DP) 


END SUBROUTINE CarbonPhenolic 
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12. Appendix E - MEM Coordinate System 


Definition of the MEM coordinate system from the MEM help file. 


COORDINATES 


All intormation displayed in the output files is relateve to the spacecratt in local spacecrat coordinates. The figure below illustrates 
the coordinate frame. 






Starboard 


A, 
-Y 


The spacecraft s velocity vector, ¥ ‘|p| , detines the X urut direction. The Y urnt direction is defined by the cross product of the 


- eyo P XP : a 
position and velocity vectors, (" “¥ IF a or the angular momentum vector, 2 The Z unit direction is then defined by the cross 
product ofthe X and Y urnt directions. 


The elevation, p, is measured from the spacecraft s velocity vector or X axis. Negative 90° points towards Earth or the nadir surface 
and positive SO? pemts tewards space or the zeruth surface, whichis racially outward along the posite 2 axis. The azimuth, 4, is 
measured positive counter clockwase trom the velocity vector or X axis through the positive Y axis or port surtace, in the XY plane. 
So O° aziniuth is the velocity vector, 90° azimuth is the port surtace or positrve Y axis and 160° azimuth is the wake surtace or 
negative acs. 


For a circular orbit, the * and ¥ vectors are perpendicular, but for elliptical orbits they are not and therefore the unit vectors for the 
body tized frame must be detined as above. For circular orbits, the mam output tle will show similar if not exact flizes and speeds 
for certain surtaces (Earth surtace and nadir, etc.). 


*e* VOTE *** The meteorcid enviranment is directional and as the spacecrajt changes directions relative to the Sun, the 


mefeoroid directionality will change aiso. Par Barth orbiting spacecrayt additional shielding by the planet will block and 
uehlock certain sources, making for rapid change in the directionality. 
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